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ABSTRACT 

The AS* = 2 6 — »• dss transition proceeds via the box - diagram in the 
Standard Model with a branching ratio calculated to be below 10~^^, thus 
providing an appropriate testing ground for physics beyond the Standard 
Model. We analyze the = 2 two - body K*-K*^, B' R-R^, 

B~ — i> K*~ and K*^ exclusive decays which are driven by the 

h dss transition, both in the Standard Model and in several extensions of 
it. The models considered are the minimal supersymmetric model with and 
without TZ parity conservation and two Higgs doublet models. All four modes 
are found to have a branching ratio of the order of 10~^^ in the Standard 
Model, while the expected branching ratio in the different extensions vary 
between lO'^ - 10"^ 



The intensive search for physics beyond the Standard Model (SM) is 
performed nowdays in various areas of particle physics. Among these, rare B 
meson decays are suggested to give good opportunities for discovering new 
physics beyond SM [|[] . Recently, it has been suggested p|, |^, H to investigate 
effects of new physics possibly arising from b — ssd ot b —* dds decays. As 
shown in Ref . [Q] , the b — >■ ssd transition is mediated in the standard model 
by the box-diagram and its calculation results in a branching ratio of nearly 
10~^^, the exact value depending on the relative unknown phase between t, 
c contributions in the box. The b — * dds branching ratio is even smaller 
by a factor of 10^, due to the relative |Vtrf/Vts| factor in the amplitudes. In 
Ref. different scenarios were used in the analysis of the b — >■ dds decay, 
which might be important in 5^ Ktt decays. The authors of Refs. [Q, Q 
have calculated the b ssd transition in various extensions of the SM. It 
appears that for certain plausible values of the parameters, this decay may 
proceed with a branching ratio of 10~^ — 10~^ in the minimal supersymmetric 
standard model (MSSM) and in two Higgs doublet models 

Thus, decays related to the b ssd transition which was calculated 
to be very rare in the Standard Model, provide a good opportunity for in- 
vestigating beyond the Standard Model physics. In Ref. [0 it was sug- 
gested that the most suitable channels to see effects of the b ssd tran- 
sition are the B~ — ^ K~K~n^ or K~ K~tt^ti^ decays. Moreover, 
when one considers supersymmetric models with 7^-parity violating cou- 
plings, it turned out that the existing bounds on the involved couplings of 
the superpotential did not provide any constraint on the b — ssd mode 
[|]. Recently, the OPAL collaboration 0] has set bounds on these cou- 
plings from the establishment of un upper limit for the B~ — > K~K~7r~^ 
decay BR{B- K-K-n+) < 1.3 x 10"^. The long distance effects in 
i?" —>■ K^K^Ti^ decay ^ have also been estimated recently and they have 
been found to be of the order 10~^^, comparable in size with the short - dis- 
tance SM contribution, thus leaving this decay "free" for the search of new 
physics. Although it appears that B~ — > K~K~tt^ or B^ —>■ K" K~tt~^7i~^ 
are very good candidates to search for the AS* = 2 transitions, we investigate 
here another possibilty for the observation of the b ssd transition: the two 
body decays of B~ . 

We consider the VV, VP, PP states. Although in principle two body 
decays would appear to be simpler to analyze, there is the complication of 

— mixing. Hence one needs also a good estimate for the b — > ssd 
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transitions as well. Nevertheless, not all the two-body states involve neutral 
K's and we shall return to this point in our summary. First, we proceed 
to describe the framework used in our analysis in which we concentrate on 
MSSM, with and without TZ parity and two Higgs doublet models as possible 
alternatives to the SM. 

The minimal supersymmetric extension of the Standard Model leads to 
the following effective Hamiltonian describing the h —>■ ssd transition ||^, 0] 

n = CMSSMis-f^dL){s-f^bL), (1) 

where we have denoted 

Cmssm = -^^[2AxMx) + QQMx)] (2) 

d 

with X = m^/mj^, and the functions fei^x) and fe^x) are given in The cou- 
plings 6fj parametrize the mixing between the down-type left-handed squarks. 
At the scale of b quark mass and by taking the existing upper limits on Sfj 
from 0] and p| the coupling Cmssm is estimated to be \Cmssm\ < 1.2 x 10~^ 
GeV~^ for an average squark mass Tn^= 500 GeV and x = 8, which leads to 
an inclusive branching ratio for b —>■ ssd of 2 x 10~^ 0. The corresponding 
factor calculated in SM is found to be 

csM = + (3) 

with f{x) and g{x, y) given in 0. Taking numerical valus from p[, neglecting 
the CKM phases, one estimates \Csm\ ~ 4 x 10"-^^ GeV'"^. 

The authors of have also investigated beyond MSSM cases by including 
R- parity violating interactions. The part of the superpotential which is 
relevant here is IV = X'^ji^LiQjdk, where i,j,k are indices for the families 
and L, Q, d are superfields for the lepton doublet, the quark doublet, and the 
down-type quark singlet, respectively. Following notations of ^ and [@] the 
tree level effective Hamiltonian is 

^ = -J2^zS^lK32K2lisRbL){sLdR) + X'^2lK32isRdL){sLbR)]. (4) 

n 

The QCD corrections were found to be important for this transition For 
our purpose it suffices to follow p retaining the leading order QCD result 
fQCD ^ 2, for = 100 GeV. 



2 



Most recently an upper bound on the specific combination of couplings en- 
tering has been obtained by OPAL from a search for the B~ — > K~K~ 7r~^ 
decay B 
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Here we take the order of magnitude, while the OPAL result is 5.9 x 10~^ 
based on a rough estimate T{B~ — >• K~K~7t~^) ~ 1/4 T{b —>■ ssd). 

The decay b — » ssd has been investigated using two Higgs doublet models 
(THDM) as well . These authors found that the charged Higgs box contri- 
bution in MSSM is negligible. On the other hand, THDM involving several 



neutral Higgses ||Tl| could have a more sizable contribution to these modes. 
The part of the effective Hamiltonian relevant in our case is the tree diagram 
exchanging the neutral Higgs bosons h (scalar) and A (pseudoscalar) 
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::^sbCsd{^ [sd) (sb) 2- (s'j^d) (5756) ) , 
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with the coupling ^j,- defined in |11] as a Yukawa coupling of the FCNC tran- 
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sitions di ^ dj. In our estimation we use the bound |^sfe^sd|/'n^ 
GeV~^, H = h,A, which was obtained in by using the Attik limit on 
C,bd/^H and assuming \C,sb/mH\ > 10^^. 



We proceed now to study the effect of Hamiltonians (|l|), (H), (||) on the 
various two body AS = 2 decays of charged B - mesons. In order to calculate 
the matrix elements of the operators appearing in the effective Hamiltonian, 



we use the factorization approximation |]T2|, |T3|, |1^] , which requires the knowl- 
edge of the matrix elements of the current operators or the density operators. 
Here we use the standard form factor representation |TB|, of the following 
matrix elements: 



{P'{p'Mrqi\P{p)) 



2 2 
2^ Tfl p Tflp, 



'-{p'-p'n, 



(7) 



where Fi and Fq contain the contribution of vector and scalar states respec- 
tively and = {p — p'Y- Also, -Fi(O) = -Fo(O) [O]. For these form factors. 
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one usually assumes pole dominance [|T3], ^ 

Fi(0) 



i^o(O) 



(8) 



and in order to simplify, we shall take my = itls- The matrix element 
between pseudoscalar and vector meson is usually decomposed as 

{Viq,ev)\qjY{l-lMPip))= (9) 



mp + my 



+ i{mp + mv 

where Q = p — q- 

2x mH + rnv 



Q 



imp + mv) 



,{p + qrA,{Q' 



2m 



V 



2my 



A2{Q') = , 



(10) 



and ^3(0) = Aq{0). For the vector and axial vector form factor we use 
again pole dominance |TB|, ITH| , and relevant parametrs are taken from III 



F(f ^(0) = 0.38, A^^*{0) = 0.32. For the calculations of the density operators 
we use the relations 



d°'{s'yab) = i{mb — ms)sh 



and 



We will use also the following decay constants: 



{Viey,q)\q,rqM=Ciq)9viq' 



and 



(11) 



(12) 



(13) 



(14) 



{Piq)\q,Yl,qm = tfpq, 

with fx = 0.162 GeV, gx* = 0.196 GeV^ Now we turn to the analysis 
of the specific modes. 
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a) B- ^ K*^K*^ decay 

For the analysis of pseudoscalar meson decay to two vector mesons it 
is convenient to use helicity formalism (see e.g. ||16[)- We denote O = 
{sj^{l — 75)0?) {sj^^l — 75)6), and then we use H = CO with C being 
1/4:Cmssm, ^/4:Csm- Using factorization and the definitions given above, 
one finds the following helicity amplitudes 

Hoo{B- ^ K*-K*'^) = CgK*{mB + mK*)[aAf''\ml,)-[3A2{ml,)] (15) 

H±±iB- ^ K*'K*') = CgK*imB + mK^)[aAf'''iml,)TlV''''''iml,)]ilQ) 
where 

with r = rriK* /ttlb, k"^ = 1 + + — 2r^ — — 2r^t^. The decay width is 
then 

T{B-^K*-K*') = JPi^[\Hoo\' + \H^+\' + \H__\']. (18) 

Within MSSM model the branching ratio becomes < 6.2 x 10~^, while SM 
gives this rate to be 6.8 x 10~^^. The TZ - parity term described by the 
effective Hamiltonian (4) cannot be seen in this decay mode when factoriza- 
tion approach is used, since the density operator matrix element {K*^\{sd)\0) 
vanishes. The two Higgs doublet model also cannot be tested in this mode 
due to the same reason. 

b) B- ^ K*-K^ decay 

The matrix element of the operator O is calculated to be 

{K%ko)K*-{k.,e)\0\B-{pB)) = -2mK*fKA^''\ml,)e* -ko (19) 
Denoting the decay amplitude by A, one finds 

EI^P = |C|VI:l^«)pA(m|,m2„mL) (20) 

pol 
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with the A(a, b, c) = + 6^ + — 2{ab + be + ac). The branching ratio is 
straightforwardly found to be BR{B~ K*' K'^)mssm < 1-6 x 10"^, which 
is comparable to the SM prediction of Ref. |T2[ for the AS" = K*~ 



decay, given as BR{B- K*- K^) = 1 x 10"^ 5 x 10"^ 2 x lO'^ obtained 
for the number of colours A^^^ = 2, A^^^ = 3, A^^^ = oo, respectively. 

The SM calculation for the AS = 2 transition leads to BR{B~ — »• 
K*-K^)sM = 1-7 X 10"^^. The MSSM which includes 7^ parity break- 
ing terms can occur in this decay. The matrix element of the operator 
Otz = + 75)c?) — 75)&) can be found to be 

{K''{ko)K*-ik.,e)\On\B-{pB)) = 

'"^^'^ -{2mK.e*-ko)A^^\ml). (21) 



{nis + md){ms + nib) 



Taking the values of the quark masses as in ||T^ nib = 4.88 GeV, = 122 
MeV, nid = 7.6 MeV and using the bound given in Eq. (|^) we obtain the 
estimation of the upper limit of the branching ratio BR{B~ —>■ K*~K^)ti to 
be 4.4 X 10~^. This limit can be raised to 1.5 x 10~^ for the upper bound on 
the couplings of 5.9 x 10~^ given in 0. 

The two Higgs doublet model gives for the amplitude of this decay 

ATHDM{B-{pB)^K%ko)K*-{k_,e) = 

o — —[2mK*fKAQ [niKje ■ ko]- 7 77 7 r, (22) 

2 ni\ (nis + nid)[nis + nib) 

which gives for the limit \S,sb^sd\/^'H > 10^^^ GeV^"^, a branching ratio of 
the order 10^^^. Due to specific combination of the products of the scalar 
(pseudoscalar) densities this is the only decay which has nonvanishing am- 
plitude within the factorization assumption. 

c) B- K-R*^ decay 

For this decay mode the matrix element of the operator O is determined 
to be 

{K*'iko,e))K-{k.)\0\B-{pB)) = 2gK*fKFf''*{nil,)e*-k^ (23) 
giving the branching ratio in MSSM with an upper limit 

BR{B- ^ K-K*')mssm = 5.9 x 10"^ (24) 
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in comparison with SM result 6.5 x 10"^'*. 

The amphtude calculated in MSSM including TZ breaking and THDM 
vanishes, due to vanishing of the matrix element of the density operator for 
K*^ state. 



d) B~ K-K° decay 

The matrix element of the operator O becomes in this case 

{K\ko)K-{k_)\0\B-{ps)) = liKFi^'imDiml-ml). (25) 

The multiplication with the corresponding 1/4Cm5sm gives the required am- 
plitude A. The branching ratio is then 

V{B-^K-K') = ^±-^Jml-Aml\A\\ (26) 

The branching ratio for MSSM is found to be BR{B^ K~ K^)mssm ^ 
2.3 X 10~^, in comparison with the 2.5 x 10~^^ found in the SM. The matrix 
element of the R parity breaking MSSM operator C*^^^ = {s'^^d){sh) is found 
to be 

{K- K'^\0^^^\B-) = {K^\s-f5d\0){K-\sb\B-) 

2 



(27) 

while the operator {s'')^h){sd) gives the same result with the opposite sign. 
The decay width is then 

T{B~^K-K'')n = — ^ v^m| - Am], \ {R-R' \0^'^\B-) /4 1 ' 



X 



f2 



i=n 



The constraint in (^ gives the bound 9.4 x 10 ^, while for the bound of 
5.9 X 10~^ for the coupling constants (P) the rate BR{B~ — > K~K^)ii can 
reach 3.3 x 10~^ 
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The long distance effects are usually suppressed in the B meson decays. 
One might wonder if they are important in decays we consider here. We 
have estimated the tree level contribution of the D{D*) which then goes 
into K{K*) via weak annihilation. We found that these contributions give 
a branching ratio of the order 10~^^ and therefore they can be safely ne- 
glected. One might think that the exchange of two intermediate states 
D{D*), K{K*) can introduce certain long distance contributions. In de- 
cay B — s> "D" " K" — i> " K" " K" the first weak vertex arises from the decay 
B " D" " K" and the second weak vertex (see e.g. 0) can be generally 
obtained from the three body decays of D ^ KKK . In Ref. it was found 
that such contributions are also very small. Therefore, we are quite confident 
to suggest that the long distance effects are not important in the two body 
AS' = 2 B decays. 



Let us turn now to the possibility of detecting these decay modes. The 
B^ — ^ K*~K*^ and B~ K~K*^ modes have clean signatures of a AS* = 2 
transition and therefore these are the channels we recommend to look for. 
The other two modes we discussed, b) and d) have a in the final states 
which complicates the possibilty of a detection because of —K^ mixing. 
Separating the desired amplitude requires the measurement of the decays of 
both Ks and Kl, since one can express [^] 



r(E- ^ K-Ks) - T{B- ^ K-Kl) 
T{B- K-Ks) + V{B- ^ K-Kl) 

where 



Re7]{B- -> KK-), (29) 



We summarize our results in the Table 1. The MSSM gives rates of the 
order 10~^ — 10~^, while the TZ parity breaking terms in the MSSM can be 
seen only in the B~ K*~K^ and B~ K~K^ decay. These are the 
modes which as we mentioned are more difficult on the experimental side. 
The THDM model can give nonvanishing contribution only in the case of 
B^ —>■ K*^K^ decay, with a rate too small to be seen. Thus, we conclude 
by stressing the possibilty of detecting physics beyond SM mainly in the 
K*~K*^, K-K*° decays. 
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Decay 


SM 


MSSM 


MSSM + 7^ 


THDM 


B- K*-K*^ 
B- K*-K^ 
B- K-K*^ 
B- K-K° 


6.9 X 10-1^ 
1.7 X 10-1^ 
6.6 X 10^1^ 
2.5 X 10-1^ 


6.2 X 10-9 
1.6 X lO-'^ 
5.9 X 10-9 

2.3 X 10-9 


10-^ - 10-6 
10-^ - 10-6 


10-11 



Table 1: The predicted branching ratio for the B~ AS — 2 two - body decays 
calculated using the factorization approach within Standard Model (the first 
column), Minimal Supersymmetric Standard Model (the second column), 
Minimal Supersymmetric Standard Model extended by TZ parity breaking 
(the third column), and Two Higgs Doublet Model (the fourth column). 
The values in columns two, three, four are upper limits, as determined from 
present knowledge of upper limits for couplings involved. 
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